AN ACCELERATED SPLITTING-UP METHOD FOR 
PARABOLIC EQUATIONS 

ISTVAN GYONGY AND NICOLAI KRYLOV 

Abstract. Wc approximate the solution u of the Cauchy problem 
d 

—u{t, x) = Lu{t, x) + f{t, x), {t, x) e (0, T] X R'^, 

u{0,x) =uoix), xeR'^ 
by splitting the equation into the system 
d 

— Vr{t,x) = LrVr{t,x) + fr{t,x), T ^\,2, ...,di, 

where L,Lr are second order differential operators, /, fr are func- 
tions of t,x, such that L = J2r-^r, f = J2rfr- Under natural 
conditions on solvability in the Sobolev spaces W^, we show that 
for any k > 1 one can approximate the solution u with an er- 
ror of order S'^ , by an appropriate combination of the solutions Vr 
along a sequence of time discretization, where 6 is proportional to 
the step size of the grid. This result is obtained by using the time 
change introduced in [7j, together with Richardson's method and a 
power series expansion of the error of splitting-up approximations 
in terms of S. 



1. Introduction 

In this paper we are interested in the rate of convergence of sphtting- 
up approximations to the solution of the parabohc, possibly degenerate, 
differential equation with time dependent coefficients 

d 

—u{t,x) = Lu{t,x) + f{t,x), G (0,T] X R"' (1.1) 

with initial condition 

u{0,x) =uoix), xgM^, (1.2) 



Key words and phrases. Cauchy problem, parabolic partial differential equations, 
splitting-up, method of alternative direction, Richardson's method. 
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where L is a differential operator of the form 

and / is a function of t > and x G M"'. The first step in the splitting 
methods is to choose suitable decompositions L = Li + L2 + ... + Ld^ 
and / = /i + /2 + ••• + /di for the operator L and the free term /, such 
that each equation 

d 

Q^Mt^ ^) = LrVrit, X) + frit, x) (1.3) 

r = 1, 2, di is integrable exactly or approximately. 

Assume for simplicity that the operators Lr and the free terms fr, 
r = 1,2, di, are time independent, and, for fixed T > and integer 
n > 1, consider the uniform step 

T^:= {t^:=tT/n,i = 0,1,2,.. .,n} (1.4) 

of step size 5 := T/n. Then a splitting-up approximation u*-"-' for the 
solution u of p.H) - p.2|) is defined by 

w„(tn := (§i'^^..§f §i'^)^Wo, ^ = 0, 1, n (1.5) 

(r) 

at the grid points. Here denotes the solution operator of equation 
fll.3p . i.e., Sj'^V is the solution of ()1.3p at time t with initial condition 
9? at t = 0. Formula ()1.5|1 means that we take ^^"^(0) = Uq, and we 
calculate the approximation at a grid point t+6 from the approximation 
u^'^\t) at the previous grid point t, by solving equations p.3|) for r = 
l,2,...,di on the same time interval [0,5] successively. First we solve 
the first equation (r = 1) on [0, 5] with initial condition Wi(0) = u^"'\t), 
and then we solve the second equation, third equation, and so on, on 
the same interval [0, 6], by taking always the value at 6 of the solution 
of the previous equation as the initial value for the following equation. 
Finally we solve the last equation (r = di) on the interval [0, 6] with 
initial condition Vd^{0) = Vdj^^i{6), and the value Vd^{6) is the value of 
the splitting-up approximation at t + 6. 

This kind of approximations is well-known in numerical analysis, 
and it has been successfully applied to various types of PDE problems. 
They are often combined with other numerical methods, as finite dif- 
ferences, finite elements, etc. Pioneering applications to the heat equa- 
tion, to hyperbolic equations, to nonlinear PDEs are presented, for 
example, in [111, El) 0, in ISH], P and in 0, [T] |2S1, respectively. 
Many applications and modifications of the splitting-up method have 
been developed in various applied fields of linear and nonlinear PDEs 
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and ODEs, under a variety of different names, like dimensional split- 
ting, operator splitting, predictor-corrector method, method of alter- 
nating directions, fractional step method, Lie-Trotter-Kato formula, 
Baker-Campbell-Hausdorff formula, Chernoff formula, split Hamilton- 
ian, split-steps, leapfrog. For guidance in the huge varieties of methods, 
names and references we refer to the survey article [12] and books [TT] , 



In the context of semigroups the splitting-up method first appears 
as Trotter's formula which can be formulated as follows: 

lim (e*^''i/"...e*^^/"e*^i/")"z = e*^z, Vz G B, 

n—*oo 

where A = Ai + A2 + ... + A^-^ and Ar are infinitesimal generators of Co- 
semigroups of contractions {e^^ '■ t > 0} and {e*^'' : t > 0} on a Banach 
space B, such that the intersection of the domains of the generators is 
dense in B. Clearly, in the context of Cauchy problems Trotter's for- 
mula states the convergence of the splitting-up approximations defined 
by the splitting 

dv 

—^ = ArVr(t), r = l,2, 

at 

to the solution of the abstract Cauchy problem 

— = Au{t), t > 0, m(0) = z. 

Our main interest in the present paper is to increase the accuracy of 
the splitting-up approximations for equation (jl.lj) . It is known that the 
error of the splitting-up approximations is proportional to 6, the step 
-size. There are, however, modifications of these approximations which 
are more accurate. A celebrated example is the Strang symmetric 
scheme 

«n(tl"^) := {^f;>,^f/rA72^f/2-^^^^^ ^ = 0, 1, n, 

whose error is proportional to 6^. This approximation scheme is pre- 
sented in ^7] , • Other symmetric schemes and their generalizations, 
are given in ^7], ^Hl? and [H]. All these schemes are of second order 
accuracy. Inspired by the above example, for given k > 2 one looks for 
a composition of splittings 

m di 

UU€s^ (1-6) 

i=i j=i 
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with real numbers c*-^ and integer m > 1 to be determined, such that 

m di 

i=i j=i 

the local error of the corresponding approximation is proportional to 
6'^'^^ in appropriate norms. Such local error leads to a global error, 
proportional to 6'', i.e., composition represents a method of (at 

least) order k. The conditions on the numbers c*-' and m which lead 
to splitting methods of high order have been studied intensively in the 
literature. Such methods are obtained in for Hamiltonian systems 
by the Baker-Campbell-Hausdorff formula. Variations of the Trotter 
formula and the Baker-Campbell-Hausdorff formula are used for linear 
and for nonlinear equations, respectively, to show the existence of meth- 
ods of any order (see [T3j, [H], [211 , [2H] and the literature therein). 
An adaptation of the method of rooted trees from the theory of Runge- 
Kutte approximations is used in ^1]. By [201 EZj? however, the 
numbers c*-^ in each scheme of order k > 3 cannot be all non- 

negative. Thus, by [201 and j2Zj the above splitting methods of order 
greater than or equal to 3 cannot be used to approximate the solution 
of partial differential equations of parabolic type. As R.I. McLachlan 
and G.R.W. Quispel write on page 392 of [T3]: "...splitting was pro- 
posed as a cheap way to retain unconditional stability. Methods with 
backward time steps can only be conditionally stable; this stumbling 
block held up the development of high-order compositions for years." 

Then the natural question arises, as to whether there exists, in the 
case of parabolic equations, a different way from the multiplicative one 
to accelerate the convergence to a higher order. One of our main results 
consists of showing that using the step size of order 6, but organizing the 
computations differently, it is indeed possible to achieve the accuracy 
of order 6'' for any k, even if Ar are (degenerate) elliptic operators 
with coefficients depending on time. In a subsequent article we intend 
to show that our method is much more universal in the sense that it 
covers very many situations in which method works and requires 
approximately the same amount of work. 

In the present paper we use linear combinations of splittings of type 
fll.Sp with different step-sizes, to achieve arbitrary high accuracy. We 
prove that for any given k >0 there exist absolute constants bo,bi, ...,bk 
expressed by simple formulas such that the accuracy of the approxima- 
tion 

Vn := boUn + biU2n + + ••• + ^fcM2'=n (1-7) 
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is of order 6^~^^ (see Theorem 12 . 21 below) . Here U2Jn is the sphtting-up 
approximation (jl.Sp along the grid (jl.4j) . with 2% in place of n. In 
particular, if /c = 1, we have to deal with two step sizes: 6 and S/2, 
and we get the order of accuracy 6"^. The Strang formula giving the 
same order of accuracy, generally, also requires working with step size 
6/2. By the way, if A = Ai + A2 and we construct our splitting- up 
scheme according to A = {l/2)Ai + A2 + (l/2)Ai, then our approxi- 
mations just coincide with the Strang one and there is no need to use 
linear combinations to get the error of order 5^. It is also worth noting 
that the above coefficients b = {b(),...,bk) are given by b := eiV~^, 
where ei := (1, 0, 0, 0) and is the inverse of the k + 1 x k + 1 
Vandermonde matrix V^^ = 2~^^^^^^^~^\ 

Our work is of purely theoretical nature and, as the referees pointed 
out, much work yet needs to be done before our results could be used in 
practical applications. We restricted ourselves to making the first step 
in attacking Problem 10 on page 492 of "For systems that evolve 
in a semigroup, such as the heat equation, develop effective methods of 
order higher than 2". However, our results show that each time when 
one has any algorithm of implementing standard splitting-up method 
to approximating the solutions of the Cauchy problem for degenerate 
parabolic equations with sufficiently smooth coefficients and free terms, 
one can improve the rate of convergence to any degree. For instance, 
we believe that usually in practice one is not doing computations with 
only one step size, and we show that having, say, three different step 
sizes each of which is of order 6 of accuracy, and just taking a linear 
combination of the results, one gets an approximation with error of 
order 6^. 

We have to admit that we do not know if our methods can be carried 
over to quasilinear equations or to equations in domains. In this connec- 
tion we note that there is a very active area of developing and applying 
in practice splitting-up methods for degenerate nonlinear convection- 
diffusion equations (see, for instance, j3] and 185 references therin). 
Our equations can be viewed as belonging to this area only if a^^ are 
constant. However, it is perhaps worth mentioning that our methods 
can be applied to solving systems of (nonlinear) ODEs and we are in 
the process of working on this subject. 

Inspired by Richardson's method we obtain our results by expanding 
the error n — n„ of the splitting-up approximation ()1.5|) in powers of 
6 = T/n. This is Theorem 12. 11 the main theorem of our paper. We use 
this expansion with 6 = T/2%, j = 0, 1, 2, k, and choose the above 
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coefficients bo,bi, ...,bk to eliminate the terms of order less than k + 1 
in the linear combination (jl.7|) . 

The main theorem of the present paper is proved by exploiting a new 
approach of [7j and [S] to splitting-up methods. As we discussed above, 
the splitting-up approximation (jl.Sp means that to get the approxima- 
tion at t + 6 from that at t, one goes back and forth in time (ii-times 
while solving equations ()1.3|) . r = 1, 2, di, successively. A basic idea 
of [7] is to arrange the splitting continuously in forward time direction, 
and to synchronize with it the original equation by time scaling. In this 
way we have differential equations for the rearranged splitting-up ap- 
proximations and for the time-scaled solution of the original equation, 
which enables us to use methods of the theory of partial differential 
equations and not semigroup theory and get an expansion for their dif- 
ference in terms of powers of 6 even if the coefficients depend on time. 
The method of [7] and [H] appeared in connection with splitting-up for 
stochastic partial differential equations. It is worth mentioning that 
most likely it is impossible to accelerate the splitting-up method in 
this more complicated situation. 

The paper is organized as follows. In the next section we introduce 
our general setting but state the results. Theorems 12.21 and 12.11 only 
for the case of time independent data for the sake of simplicity of 
presentation. Theorem 12.11 is proved immediately after its formulation 
on the basis of Theorem 12.21 which in turn is proved in Sec. |31 after we 
prepare some auxiliary facts in Sec. El In SecElwe generalize Theorems 
12. 21 and im for time dependent data and derive some consequences valid 
in the time-homogeneous case as well. 

In conclusion we introduce some notation used everywhere below. 
Throughout the paper d > l,di > 2 are fixed positive integers, K,T 
are fixed finite positive constants, and 

A := d/dx\ Dij ■= d'^/dx'dx\ Dt := d/dt. 

We denote by the Sobolev space defined as the closure of 
functions (/? : M'^ M in the norm 




where D"* := Dj^ ...Dj'' for multi-indices 7 = (71, ...,7d) of length I7I : = 
7i + 72 + ••• + 7d- Unless otherwise indicated, we use the summation 
convention with respect to repeated indices. 

We express our sincere gratitude to the referees for useful criticism 
which helped improve the presentation. 
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2. Formulation of the main results. The case of time 
independent coefficients 

We consider the problem 

Dtu{t,x) = Lu{t,x) + f{t,x), t e (0,T], X e M'^, (2.1) 

u{0,x) =uo{x), xeR'^, (2.2) 

where L is an operator of the form 

L = a'^it, x)Dij + a\t, x)Di + a(t, x), 

/ and Mo are real functions of (t, x) G (0, T] x and of x G M"^, respec- 
tively. We assume that the coefficients a*-', a*, a and the derivatives a^^ 
of a*-' are bounded Borel functions of (t, x). We fix p > 2 and assume 
that Mo and / are measurable and |mo|^ and |/|^ are integrable over 
and over [0,T] X R"^, respectively. 

Definition 2.1. By a solution of problem ()2.H) - ()2.2|) we mean an W^- 
valued weakly continuous function u{t) = u{t, ■) defined on [0,T] such 
that for all G CoiR'^) and t G [0, T] 

(M(t, ■),</)) = (u(0, ■), 0) + / [-(a*^' An(s), Z},0) 

Jo 

+ ((a* - a^j)Diu{s) + au(s) + /^(s), </>)] (is, 

where ( , ) denotes the usual inner product in L^(]R'^). Quite often we 
write equation (j2.H) and similar equations in the form 

du{t) = {Lu{t) + fit)) dt 

bearing in mind the differential of u in t only. 

Suppose that we split equation ()2.1|) into the equations 

Dtv{t,x) = Lrv{t,x) + fr{t,x), t G (0, T] , X G (2.3) 

with 

di di 

Lr := al^{t,x)Dij + al{t,x)Di + ar{t,x), L = ^Lr, f = ^fr, 

r=l r=l 

such that these equations are more pleasant from point of view of nu- 
merical methods than the original one. This motivates the multi-stage 
splitting method, which we describe below. First we need some as- 
sumptions. 

Fix an integer / > 1. 
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Assumption 2.1 (ellipticity of L^). For each r = 1, 2, for dtxdx- 
almost every (t, x) G [0, T] x 

al^{t,x)yX^ > 

for all (A\A2, ....A"^) G R'^. 

Assumption 2.2. (i) The partial derivatives 

D^DPai^, D'tDPai, D^DPar for z, j = 1, 2, ci, r = 1, 2, di 

exist and by magnitude are bounded by K for all integers s > and 
multi-indices p, satisfying 2s + \p\ < I. 

(ii) For every integer s G [0, //2] 

sup \\Dtfrmi-2s,p<K. 

te[o,T] 

(iii) We have Uq G and < K. 

It is well-known that under the above conditions equations (|2.H) 
and (|2.3|) with initial condition m(0) = Uq admit a unique generalized 
solution u and v, respectively, which are PV^- valued weakly continuous 
functions of t > (see, for instance. Theorem 13. II below). We want to 
approximate the solution u, by using the splitting- up method, i.e., by 
solving equations ()2.3|) successively with appropriate initial conditions 
on appropriate time intervals. Let us formulate now our splitting-up 
scheme in the case when the coefficients a*-', a* , and free terms fr are 
independent of the time variable t. 

Set Tn := {ti := iT/n : i = 0, 1, 2, n}, 6 := T/n for an integer 
n > 1. Then for fixed n we approximate the solution u of ()2.1|) - ()2.2j) 
at ti = iT/n recursively by Un{0) := uq, 

Un{ti+i) := Sf'l.Bf^S's^UniU), ^ = 0, 1, 2, n - 1, (2.4) 

where S|^V •= '^(^) denotes the solution of equation ()2.3p for t > 
with initial condition v{0) = tp. 

It is known that if Assumptions 12. l[ 12.21 are satisfied with I = 771 + 4, 
then 

max \\u(t) - Unit)\\m,p < N/n 

for all 77, > 1, where N depends only on d, di, T, K, p, m. Moreover, this 
rate of convergence is sharp (see |H| , where this result is a special case of 
the rate of convergence estimates for stochastic PDEs). In the present 
paper we want to show that by suitable combinations of splitting-up 
approximations we can achieve as fast convergence as we wish. We 
show this by the aid of the following theorem on expansion of m„ in 
powers of the step-size 5. 
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Theorem 2.1. Let m > and k >0 be integers. Let Assvmptions \2. 1\ 
and \^.^ hold with 

l>A + m + Ak. (2.5) 

Suppose that the coefficients and the free terms fr do not 

depend on t. Then for all n > 1 and t ^ Tn and x G W^, the following 
representation holds 

Un{t, x) = u{t, x) + Su^^\t, x) 

+6'u^'\t,x) + ... + 6'u^'\t,x)+R'i\t,x), (2.6) 

where the functions u'^^\ . . . ,u^^\ and R^n \ defined on [0,T], are W^- 
valued and weakly continuous. Furthermore, u^^\ j = l,2,...,k, are 
independent of n, and 

sup||i?W(t)|U,,<iV5'^+^ (2.7) 

for all n, where N depends only on I, k, d, di, K, m,p, T. 

Remark 2.1. If k = and p = 2, the resuh holds under a weaker 
restriction on /: / > 3 + m (see, for instance j7]). For general p > 2 
and k = the result is proved in [Sj. 

We prove Theorem 12.11 in Section HJ Now we deduce from it a result 
on the acceleration of the splitting-up method. Let V denote the square 
matrix defined by V'^ := 2~'-'-^^^^-^\ i,j = l,...,k + 1. Notice that 
the determinant of V is the Vandermonde determinant, generated by 
1,2^^, ...,2^'', and hence it is different from 0. Thus V is invertible. 
Set b := {bo, 6i, bk) := (1, 0, 0, 0)V-\ and define 

k 

Vn(t) := ^ bjU2Jn(t), t ^ := {iT/n : i = 0, 1, .., n}, 

j=0 

where U2jn is the splitting-up approximation based on the grid T2i„ : = 
{iT/(2%) : i = 0,1, ..,2%}. 

Theorem 2.2. Let m > and k > be any integers. Let Assumptions 
\2.1\ and \2.^ hold with I satisfying Suppose that the coefficients 

al? , a* , ttr and the free terms fr do not depend on t. Then 

max||t;„(t)-n(t)|U,p<iV5'=+\ 
where N is a constant, depending only on I, k^d^di^ K^m^p^T . 
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Proof. By Theorem 12.11 

k 



U2.n = u + Y^ — n« + J = 0, 1, k. 

i=l 

Therefore for all > 1 

A: 
j=0 



U 



since X]j=o ~ ^^"^ Sj=o = for z = 1, 2, ...fc by the definition 

of (6o, hk). Hence = ^^^^^ ^j^a^i' ^'^^ 

max||^;„(t) = max || ^ 6ji?gJ^(t) 

" i=o 



bjU2in — 


k 

i=o 




••+1:''- 

i=o 


1=1 


i=o 




i=o 



<^|6,|max||i?gi(t)|U,,<iV5^+\ 

j=0 

by ()2.7|) . where is a constant depending only on /, T, K., d, di,m,p,k. 

□ 

Remark 2.2. Assume that u'-^-' = in expansion ()2.(ij) . This hap- 
pens, for example, for Strang's splitting, which is a special case of 
our splitting-up scheme, as it is explained in the Introduction. In this 
case we need only take k terms in the linear combination to achieve 
accuracy of order k + 1. Namely, we define now Vn,{t) by 

'"nit) ■■= \jU2Jn{t), t E Tn, 

j=0 

where 

(Ao,Ai,...,Afc_i) :=(l,0,...,0)\/-\ 

and V is now a. k x k Vandermond matrix with entries Vn := 1, Vij : = 
2-(*-i)j for i = l,2,...,/c and j = 2,...,k. Then Theorem 12.21 remains 
valid, what one can prove in the same way as Theorem 12.21 is proved. 
For example, 

1 4 

Vn{t) := --Un{t) + -M2„,(t), t G 

is an approximation of accuracy 6^ in the case of Strang's splitting. 
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3. Auxiliary Results 
Let us consider the partial differential equation 
du{t,x) = {Lu{t,x) + f{t,x))dA{t), te(0,T], xeM^ (3.1) 

u{0,x) =uo{x), xeR'^, (3.2) 
where L is an operator of the form 

L = a'^it, x)Dij + a\t, x)Di + a(t, x), 

A = A{t) is a continuous increasing function starting from 0, / and Uq 
are real functions of [t, x) G (0, T] x R*^ and of a; G W^, respectively. Fix 
an integer I > and a real number p > 2. We understand the solution 
in the spirit of Definition 12.11 and make the following assumptions. 

Assumption 3.1 (smoothness of the coefficients). The coefficients of 
L are measurable. The derivatives in x G M"^ of the coefficients a^^ up 
to order 2 V /, of the coefficients a*(t, x) up to order 1 V Z, and of a{t, x) 
up to order / exist for any t G (0, oo), and by magnitude are bounded 
hjK. 

Assumption 3.2. We have 

uoeW^, f eLp{[0,T],wl). 

Assumption 3.3 (ellipticity of L). For alH > 0, x G M"*, and A G R'^, 
we have 

a'^{t,x)X'\^ > 0. 
Assumption 3.4. The function A is absolutely continuous and 

A{t) := jA{t) < K 
for (it- almost every t > 0. 

The following result is well-known in PDE theory (after replacing 
dA in p.lj) with Adt we easily get it, for instance, from or from 
Theorem 3.1 in [Zj). 

Theorem 3.1. Under Assumptions \3. 11 [Q. [Q. and\3.4\ with I > 1 
the Cauchy problem ^^^-^^^ has a unique generalized solution u. 
If Assumptions Vj.lV VJ.'A \3.'J[. and \3.4\ hold with I > 0, and u is a 
generalized solution of ^^H})-^^^, then for every integer li G [0, /] 

sup Mt)\\l^<N{\\u4l^+ f\\fmi,dt 

where N is a constant depending only on T, K,l,p,d. 
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Under the assumptions of Theorem 13. II let TZf denote the solution of 
equation (j3.1|) with initial data uq = 0. Then by virtue of Theorem 13. II 

'}Z:Lp{[0,T],Wl)^C^{[0,T],Wl) 

is a bounded linear operator, where C^([0, T], Wp) denotes the Banach 
space of weakly continuous M^p- valued functions u = u(t), t G [0,T] 
with the norm sup^gjQ ^-] ||M(t)||i^p. 

Let us now consider the equation 

du{t, x) = Lu{t, x) dA{t) + g{t, x) dH{t), (3.3) 

{t,x) e (0,T] X M'^, 

where g is a. real-valued function of (t,x) G [0,T] x M.'^ and H is an 
absolutely continuous function of t G [0,T]. 

Assumption 3.5. We have g G Lp([0,T], Wj,'^'^), and there exists g' G 
Lp{[0,T],Wl) such that 

d{g{t),(j)) = {g'{t),(j))dA{t), tG[0,T] 

for all G Co(M'^). 

Lemma 3.2. Under Assumptions \y. 11 [Q. \3.4\ [Q with I > 1 egna- 
tion !^} with zero initial data has a unique generalized solution u. 
Moreover, 

u = n{H{Lg - g')) + Hg =: Q{H, g). (3.4) 

// Assumptions IJ.il [Q. \3.4\ [Q hold with I > and equation \3. 'J\) 
with zero initial data admits a generalized solution u, then for every 
integer li G [0, /] 

sup \\u{t)\\i^^p<N sup \H(t)\[ sup \\g{t)\\i^^p 

te[0,T] te[0,T] ^tG[0,T] 

+ {j\hmU2,r>+\\9'mi,,)dtY''). (3.5) 

where N is a constant depending only on p,d, K,l,T . 

Proof. Note that 

dMt), ct))H{t)) = {git), 0) dHit) + ((7'(t), dAit), 

for all G Co(M'^). Therefore u solves equation ()3.3|) . when U7(t, x) : = 
u{t, x) — H(t)g(t, x) solves 

dw{t,x) := {Lw{t,x) + H{t){Lg{t,x) - g'{t, x))} dA{t). 

Hence equality ()3.4|) follows by Theorem 13.11 and it implies ()3.5|) . □ 
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Remark 3.1. We often consider equation (j3.3|l when is bounded 

Then under Assumptions 13.11 13.31 13.41 13.51 with / > 1, equation (j3.3|) 
with zero initial data has a unique generahzed solution u by Theorem 
13.11 By Theorem 13.11 this solution belongs to C^([0, T], VTp) and its 
norm in this space admits an estimate with a constant depending on 
the bound for H and only the lvp([0, T], W^-Jioim. of g. It is important 
that the function H enters ()3.5p only through sup \ H\ and not any char- 
acteristic of its derivative, however for that we pay a price requiring g 
to have more derivatives. 

4. Proof of Theorem 12.11 

Throughout this section the assumptions of Theorem 12.11 are sup- 
posed to be satisfied. In particular, / > 4. Fix n and introduce 5 = T/n. 
We use the idea from [7] and |H] of rearranging the splitting method in 
forward time. We achieve this by considering the equation 

di 

dw{t, x) = '^^{Lrw{t, x) + fr) dAr{t), w {0, x) = Uo{x) , (4.1) 

where the time change A^, r = 1, ...,di, is defined by the requirements 
that Ar{0) = 0, Ar(t) be absolutely continuous, and its derivative in 
time Ar be periodic with period di6 and 

= l[,_i,,](t/5), te[0,di5] (a.e.). (4.2) 

Instead of the original Cauchy problem ()2.1|) - ()2.2j) we consider 

dv(t, x) = {Lv(t, x) + /) dAQ{t), f (0, x) = Uo{x), 

where 

Aoit) := t/di. (4.3) 
Clearly, v{t) = u{Ao(t)), and 

v{dit) = u(t), w{dit) = Unit) for all t e T„. 

Therefore our aim is to show that Theorem 12 . 11 holds with v and w in 
place of u, and Un, respectively, for all t = idi6, i = 0, 1, ...n. To this 
end first we introduce some notation. We call a sequence of numbers 
a = aia2.--0ii a multi-number of length |a| := i, if aj G {0, 1, 2, di\. 
The reader should notice the difference between multi-numbers and 
multi-indices. The set of all multi-numbers is denoted by M. For every 
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multi-number a we define a function '■ [0, oo) R and a number 
Cq, recursively starting as follows: 

5^ :=(5-^(A^-Ao), c^ = for 7 = 0,l,2,..rfi. (4.4) 

If for every multi- number /3 = (3i...(3i of length i the function 5^ and 
the number cp are defined, then 

cp^ ■= 5"^ / Bf3{s)A^{s) ds, (4.5) 
5^^(t) := [ {Bp{s)A^{s) - c^^io(s)) ds (4.6) 

JO 

for 7 = 0, 1, 2, (ii, where A^{s) := (i/l^(s)/(is. 
Notice that by ()4.6|) we have 

(t) dA^ (t) = Cf3^ dAo {t) + 6 dBp^ (t) (4. 7) 

for all multi-numbers (3 and 7 = 0,1,2, ...,di. We will often make use 
of this equality and of the following lemma. 

Lemma 4.1. For every a G A/" the function B^ is di6-periodic, i.e., 
Ba{t + di6) = Ba{t) for all t > 0, and Ba{idi5) = for all integer 
i > 0. Moreover, the numbers c^, the functions Ca{t) := Ba{St), and 

sup \Ba{t)\ = sup \Ca{t)\ 
t>0 t>0 

are finite and do not depend on 6. 

Proof. That the first assertion is true for a = 0, is almost obvious. 
If it is true for a = f5, where is a multi-number, then the integrand 
in ()4.6p is (ii(5-periodic and by definition of Cf^^ its integral over the 
period is zero. It follows that the first assertion holds for a = /S'j, so 
the induction on the length |a| finishes the proof of the first assertion. 

To prove the second one we again use the induction on the length 
i = \a\. This statement is true when |a| = 1. Assume that it is true 
for all multi-numbers f5 of length i and notice that according to (j4.2j) 
and ()4.3|) A^{6s) are di-periodic in s and independent of 6. Therefore, 

Cf3^ = - / Bi3{s)A.y{s) ds,= / Ci3{s)A^{6s) ds 
^ Jo Jo 

is independent of 6 by the induction hypothesis. Similar argument 
works for Ca{t). □ 
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We use the notation TZf and Qag, oi G A/", for the solutions of equa- 
tions (j3.1|) and (j3.3|) . respectively, with zero initial condition and Aq 
and Ba in place of A and if, respectively. Notice that, unlike in the case 
of uniformly parabolic operators, TZ and do not increase regularity. 

The following lemma exhibits our two main technical tools: centering 
Ba and integrating by parts with respect to t. 

Lemma 4.2. Take some functions 

heLp{[Q,TlWl), /i, GLp([0,r],W;i), r = l,...,rfi, /io = 0. 

Let u he a solution of the "equation" 

di 

du = hr dAr 

r=l 

m(0) G W^, which is a particular case of equation when Lr = 0. 

Finally, let Lu G i^p([0, T], W^). Then for any a E N 

TZ{B^h) = Caonh + 6Qo,oh, (4.8) 

QaU = Tl{CaoLu — CarK) + SQaoLu — SQarK + B^U. (4.9) 

Proof. To prove ()4.8|1 it suffices to use the definitions of TZ and Q/3 
(see Theorem 13.11 and Lemma I3.2|l and use that by virtue of ()4.7j) for 
ip = 7l{Bah) we have 

dip = Lip dAo + Bah dAQ = Lip dAo + Caoh dA^ + 6h dB^Q. 

To prove ()4.9p observe that by definition 9 := QaU satisfies 

d9 = LOdAo + udBa, ^(0) = 0. 

This and (j4.7p imply that ip := 6 — uB^ satisfies V'(O) = and 

dip = Lip dAQ + LuBa dAQ — h^.Ba dA^ 

= Lip dAo + C^qLu dAQ + 5Lu dB^O — Carhr dAQ — 6hr dBar- 

Now ()4.9p follows from the definitions of TZ and Qp. The proof of the 
lemma is complete. □ 

Now we introduce some differential operators L^ and functions 
defined for multi-numbers 7 as follows: Lq := 0, /o := 0, 

Lj . Lr, f^ . /j. 

for 7 = r G {1, 2, di}, and 

L^Q \— LL^, L^j. := LjL^ 
/7O • Lf^j f^r ■ L^fr 
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for r = 1,2, ...,di. 

In this notation we have the foUowing. 

Lemma 4.3. Let a, p E M and 2\p\ + 3 < I. Then 

QaiLpW + fp) = Car'R-iLfSrW + ffSr) 

+ 6QariLprW + f^r) + B^{LpW + fp). (4.10) 

Proof. It follows from formula ()4.9|) applied to u := Lpw + fp, when 
hr = LjsLrW + Ljjfr (remember fr are independent of t), that the left 
part of ()4.1U|) equals 

which is easily seen to be equal to the right-hand side of ()4.1Up . 

□ 

We derive from (|4.1(Jj) one of the most important formulas. 
Proposition 4.4. Let n > be an integer and / > 2k + 3. Then 

K 

W = V + ^6'Y^ Ba{L^W + fa) 
i=l \a\=i 

K 

+ E c^niLaW + fa) + S''+V^'> (4.11) 

i=l \a\=i+l 

for all t G [0, diT], where 

q|=k+1 

Proof. First notice that foT (po := w — v we have 

dipo = d{w — v) = LipQ dAo + 6{LrW + fr) dBr, 

which proves (|4.1ip for k = 0. 

Next we fix a K > 1 and transform r^'^\ for i = 0,...,K— l,by applying 
(ICTlll with a = /3 and |a| = i + 1 when e W^p^^'+l Then we get 

r^'''' = ^ Cap7l{Laf3W + /a/3) +5 ^ Q{LapW + /^/j) 

|a|=j+l,|/3|=l |a|=i+l,|/3|=l 

+ ^ Ba{LaW + fa)= ^ Ba{LaW + fa) 
\a\=i+l |a|=j+l 

+ Yl CaUiLaW + fa) + 6r^'+'\ 
\a\=i+2 
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This shows how r^^\ r^^\ r^'^^ are related to each other and certainly 
proves the proposition. 

□ 

Decomposition ()4.1H) looks very much like ()2.6|) the only difference 
being that the factors of 6^ depend on the approximating function w 
and the coefficients of 5^ in the second term on the right contain Ba 
which is no power series in 6. However, observe that we have to estimate 
the difference v — w only at the points idi6 at which all vanish. 

Our next step is to "solve" (|4.1ip with respect to w by the method 
of successive iterations, that is by substituting w given by (j4.11|) into 
the right-hand side of the same equation. In the process of doing so 
we encounter only one difficulty when the second term on the right 
is plugged into the third one and we have to develop expressions like 
lZ{BaU) into power series in 5. We transform these terms by using ()4.8|) 
and 

First we introduce the notation 

wp = Lfsw + ffs, 
observe that in these terms ()4.10|) is rewritten as 

QaWf3 = CarTlWpr + B^^Wp + 5QarW^/3r, (4-12) 

and note the following. 

Lemma 4.5. If n > is an integer and a,(3 E M and 2(|/3| + fi;) + l < /, 
then 

K 

TZ{BaWp) = X] CaO-yTZWi3^ 

i=0 |7|=i 

+ BaO-yWf3^ + 6^+^ J2 Q^O-yWfS-y, (4.13) 

i=l 171=4— 1 |7|=K 

where for any multi-numbers fi, v 

Cy^TZvOf,^ := CuTZw^, ^ B^^^w^^ := B^^w^, 

|7|=0 |7|=0 
|7|=0 

Proof If K = 0, ()4.13p follows from (fl^ . If k > 1, by applying 
repeatedly ()4.12j) as in the proof of Proposition 14.41 we find 

ft-i 

i=0 |7|=i+l 
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i=0 |7|=i |7|=K 

We use this formula for aO in place of a and finish the proof by referring 
to 

□ 

Let M. denote the set of multi-numbers 7i72--.7i with 7j G {1,2,..., di}, 
i = 1,2, and integers i > 1. 

Lemma 4.6. The following statements hold. 

(i) Let 7 = 7i72--.7i & M. he such that |7| = z < 1 + 1/2. Then 

L-y = (—1)''''' ^L^^...L^-, = (—1)''^' ^ L^^...L^._^f^., 

(ii) Let e M he such that \(3\ + I7I < 1 + 1/2. Then 

LpLj = —Lpj, Lpf^ = —fp^y. 

(in) Let a E M he such that p := |a| < 1 + 1/2. Then there exist 
constants 0(7) = c(a,7) G {0, ±1} defined for all E with I7I = p, 
such that 

La= Yl f»= E (4-14) 

7e7V1,|7|=p 'y€M,\j\=p 

Proof. Part (i) follows immediately from the definition of L^, f^ by 
induction on |7|. Part (i) obviously implies Part (ii). Part (iii) clearly 
holds for a = and a = r G {1, ...,di}. Assume that equations ()4.14p 
hold for some a G Af, |a| < 1 + 1/2. Then 

Lar = -LaLr = - ^ c{'j)L^Lr = ^ c{'j)L^r, 

7| = |a| j7| = |Q! 

far = -Lafr = " J]] c{-f)L^fr = ^ 0(7)/^^ 

7| = |a| |7| = |q^ 

for r G {1, 2..., di}, and 

di d\ 

r=l 7gX,|7|=p r=l 7gX,|7|=p 

d\ d\ 
faO = Lfa = 'YLr ^ c{-f) f^ = - ^ ^ c{-f)fr^, 
r=l 7eX,|7|=p r-1 ■y^M,\j\=p 

which prove (iii) by induction on |q!|. □ 
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We introduce sequences a = {Pi, P2, A) of multi-numbers Pj G 
where z > 1 is any integer, and set \a\ := \Pi\ + \P2\ + ••• + We 
consider also the 'empty sequence' e of length |e| = 0, and denote the 
set of all these sequences by J'. For a = {Pi, P2, Pi), i > 1, "we define 

Sa = TlLp^ ■ ... ■ TZLp^ 

and for a = e we set 

Se = n. 

Notice that involves 2\a\ derivatives with respect to x and certain 
number of operators TZ which do not increase regularity. Therefore, 
basically, S^j has the power of a differential operator of 2|cr|'th order. 
If we have a collection of functions g^, indexed by a parameter u taking 
values in a set A, then we use the notation 

for any linear combination of Qi, with coefficients independent of the 
argument of and of 6. For instance, 

A (o-,7)eA (o-,7)eA 

where c{a, 7) are certain constants independent of 6. These constants 
are allowed to change from one occurrence to another. 

For functions u = u{t,x) = u{6,t,x) depending on the parameter 6 
we write u = Om{S'^), if 

sup^"'' sup \\u{t)\\m,p < 00. 

S te[0,diT] 

We also use the following sets 

A{{} = {{a,P) ■.aeJ,PeM,\a\ + \P\ < z}, 

B{i,j) = {{a,P) ■.aE^^,PeM,\a\=^, \P\ < j}. 

Lemma 4.7. Let n, fi > be integers and a G A/", /3 G M., o E J . 
Assume that 

2(|(t| + |/?| + +/i + 2 < /. (4.15) 

Then 

K ^ 

i=0 A{\a\ + \^\+i) 

K 

+ E^^ B^,wp,+0,{6^+'). (4.16) 

1=1 B{\a\+i,\a\ + \l3\+i~l) 
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Proof. For a = e, when = TZ, equation (j4.16p turns out to be just 
a different form of (j4.13p . wliicli is applicable since 2{\P\ + k) + 1 < I. 
Indeed, owing to Lemma 14.61 (iii) 

\l\=i MH + W+i) 

|7|=j-l B{\a\+i,\a\ + \l3\+i-l) 

Furtliermore, for |7| = k (see Remark 13. ip 

Qao-yWp^ = 0^{1), since 2{\f]\ + n) + fi + 2 < I. 

For \a\ > 1 we proceed by induction on the length ^(5*0-) of 5*0- = 
TZLp-^ ■ ... ■ TZLp., which we define to be j. If i{Sa) = 1, then 5*0- = TZL^, 
for a z/ G 7W with u = a and it suffices to notice that 

S^iB^wp) = nL^^BaWfs) = -TZ{BaWyp) = -Se{B^wpi), (4.17) 

where (3' = v(3 e M a nd 2(|/3^| +Av) +/x + 2 = 2(|a| + +^ + 2 < I. 

Assume that ()4.16|) holds whenever l{So) = s, and take an such 
that £{Sa) = s + 1. Then 5*0- = TZLi^Sa', where z^, cr' G Ai, |i^| + |(t'| = \a\ 
and i{Sa-') = s. Furthermore, 

2{\a'\ + |/?| + +/i' + 2 < /, 

where n' = /x + 2|z/|. By the induction hypothesis 

K ^ 

Sa'iBaWp) = '^5' ^ Sa^Wf3^ 
i=0 A{\a'\ + \l3\+i) 

K 

+ E* B^,wp,+OAS''^'). 

i=l B(|Q|+i,k'| + |/3|+i-l) 

We apply TZL^, to both parts of this equality and take into account that 
LyWp^ = —Wyp^ and + \a'\ = \a\. Then similarly to ()4.17|) we get 
that 

SaiBaWfs) = ^S' ^ S^^WfS^ 
i=0 A{\a\ + \f3\+i) 

K 

+ 5Z* S,{B^,wp,) + 0,{5^^^). (4.18) 

i=l B(|a|+i,|o-| + |/3|+j-l) 

Now we transform the second term on the right. Take (ai,/?i) G 
B{\a\+i,\a\ + \l3\+i-l) and notice that then < |(t| + + « - 1. 
Hence by assumption ()4.15p 

2{\(3i\ + K-i)+ii + 2< I. 
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Therefore, by the result for cr = e 

K—i ^ 

j=0 A{\M+j) 

K—i 

j=l B{\ai\+j,\(3i\+j~l) 

We substitute this result into ()4.18|1 and obtain ()4.16|1 after collecting 
the coefficients of and noticing that, if [ai, /3i) G + i, jo"! + 

\(3\ +1-1) and (a2,/52) G B{\ai\ + j, \/3i\ +j- 1), then 

(a2, (32) eB{\a\+i+ j, |a| + + z + j - 1). 

This justifies the induction and finishes the proof of the lemma. □ 

We remind the reader that throughout this section the assumptions 
of Theorem 12. II are supposed to be satisfied and in the following propo- 
sition use the notation 

i 

B\i,3)= \jB{ir,3). 
n=i 

Proposition 4.8. For any j = 0, 1, k we have 

j k 

w = v + j2^'J2*^-^f^ + 5Z J2* -^-i^/^i 

i=l A{2i) i=j+l A{i+j+l) 

k 

+ J2* B^.wp.+OUS"^'), (4.19) 

i=l B*(i,i+j) 

where := Lpv + fp. 

Proof. By Proposition 14.41 (notice that, due to ()2.5|) . / > 2A; + 3 and 
2(A; + 1) + m + 2 < /) we have 

k k 

i=l \l3\=i i=l \l3\=i+l 

which means that ()4.19j) holds for j = 0, since by Lemma [4.61 (iii) 

|/3|=i |/3|=i 7e7H,|7|=j B*{i,i) 

|/3|=i+l l/3|=*+l 7GX,|7l=*+l A(i+1) 
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Next, assume that k > 1 and (j4.19|) holds for a j G {0, k — 1}. 
Transform the first term with z = j + 1 in the second sum on the right 
in ()4.19|) by using Lemma 14.71 To prepare the transformation take 
(cri,/3i) G A{2i) = A{i + j + 1) so that \ai\ + \f3i\ < 2i and apply the 
operator S^^L/^-^ to both parts of equation ()4.1H) with k — i in place of 
K. Then we obtain 

k—i 

ii=l |ai|=ii 

k—i 

n=l |ai|=ji+l 

where 

\a\=k-i+l 

Owing to 

l-2{k-i + l + \(3i\ + \ai\)>l-2{k + i + 1) 
> l-2{2k + l) >m + 2, 

we have r^'^"*) = Om(l)- By the way, this is the only place where we 
need I to be not smaller than 4 + m + Ak. Hence by Lemma (4.61 (iii) 

k—i ^ 

n=i (a2,/32)e-B(n,|/3i|+n) 

k—i 

+ S^,wp,+0m{5''-'^^)=:J^ + ... + J,. (4.20) 

n=l A(|(Ti|+|/3i|+ji+l) 

Now Lemma 14.71 with k — i — ii in place of k and m in place of /i 
allows us to transform terms entering J2. For |/32| < lAI + ^1 we have 
(remember that ((Ti,/3i) G A{2i)) 

2(|ai| + \p2\ + k-i-ti) + m. + 2< 2(|ai| + + k - i) + m + 2 

< 2{i + k) + m + 2 < Ak + m + 2 < I. 

Therefore 

k—i—ii ^ 
S„^{Ba^Wf3^) = Y 5Z 

42=0 A(|cri| + |/32|+j2) 

k—i—ii 

+ E E* 5.3^^+0^(5'=-^-^^+^). 

12 = 1 B(|a2|+i2,ki| + |/32|+i2-l) 
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We plug this result into J2 and in order to collect the coefficients of 
§ii+i2 notice that, for (a3,/53) e A{\ai\ + [/Jsl + ii) and (a2,/32) G 
B{ii, \[3i \ + ii) it holds that 

kal + I/53I < kil + I/52I + i2 < kil + + ^1 + ^2- 
Furthermore, if {a^,l3^) G -B(|a2| + ^2, + 1/^21 + «2 — 1), then 

ksl = |a2|+^2 = ^1 + ^2, l/^sl < 10-11 + |/52|+«2-l < |0'l| + |/5l|+«l + «2- 

It follows that J2 is written as 

k—i 

ii=l A(|cri| + |/3i|+ii) B(h,\<Ti\ + \l3i\+n) 

which just amounts to saying that visually in the definition of J2 one 
can erase S^^-, carry all differentiations in it onto wp^, and still pre- 
serve fionii . 

Then we see that 

A(2j+2) A{2j+2) 
n=l A(|<Ti| + |/3i|+ji+l) B(ii,|(7i| + |/3i|+ii) 

Next we notice again that, for ((Ti,/5i) G A(2j + 2) and |cr2| + |/52| < 
Ic^il + |/5i| + zi + 1, we have |cr2| + |/32| < j + 2 + zi + j + 1, whereas if 
\h\ < lo^il + + «i, then I/32I < j + 1 + «i + j + 1- Therefore, after 
changing ii + j + 1 i (> j + 2) we get 

A(2j+2) A(2i+2) 

+ E '^'( E '^'^2^/32+ E B^^wp^). 

i=j+2 A{i+j+2) B"{i,i+j+l) 

This shows that the term with z = j + 1 in the second sum on the right 
in ()4.19p can be eliminated on the account of changing other terms 
with simultaneous shift j ^ j + 1. Thus the induction on j proves the 
proposition indeed. □ 

Now we finish the proof of Theorem 12.11 By taking j = in Propo- 
sition UHl we find 

k 

yj = v + J2s'w^'^ + E 4a,p,S)B^wp + 0m{S''+^), (4.21) 

4=1 B*{k,2k) 
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where 

:= J2*SaVf3 e a([0,T], W^), t = l,2,...,fc, 

A{2i) 

are independent of 6, and c(a, /3, 5) are certain constants. It is not hard 
to follow our computations in order to see that 

sup sup r('^+i)||0™(5'=+i)(t,-)IU,p<iv, 

te[0,diT] n,S=T/n 

where the constant depends only on d, di, T, K, k, m, p, I. After that 
to finish the proof it only remains to recall that Ba{jdiS) = for all 
integers j > and 

v{dit) = u(t), w{dit) = Unit) \/t E Tn = {iT/n : i = 0,1,2, ...,n}. 

5. The case of time dependent coefficients 

We consider here the Cauchy problem (jl.ip - (|1.2j) for time dependent 
coefficients. We split, as before, the coefficients and the free terms into 
di terms, 

di d\ d\ di 

r=l r=l r=l r=l 

define 6 = T/n, ti = = 5i, = {U : i = 0, 1, n}, and keep 
Assumpt ions l2.1ll2.2[ As before, we also denote Lj. := a^J Dij+a^Di+ar- 
One of splitting-up approximations Unit) for t G T„ is defined as 
follows. Let S^j'' be the operator mapping each function Lp of an appro- 
priate class into the solution of the problem 

Dtvit,x) = diLrvit,x) + difr{t,x), t > s, vis, x) = (fix). 

Then the approximations are introduced according to 

M„(0) := Uq, 



Unit,+i) := ^Xr-^LuXL^niti), 2 = 0, 1 , 2, . . . , n, (5.1) 

where tij := ti + jS/di, for j = 1, 2, (ii — 1, J := (ii — 1. 

There are many other ways to extend the splitting-up approxima- 
tions ()2.4j) to PDEs with time dependent data. Along with ()5.H) we 
also consider another approximation, which has the advantage that in 
each step we need to solve a time independent PDE, which is usually 
more convenient in practice than solving time dependent PDEs. This 
time we define the approximation u„ by 

M„(0) := Mo, 
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«„(tr+i):=Si''^(^m)-Sf (^m)Si'^(^mK(tr), ^ = 0, 1, 2, n, 

(5.2) 

where E'g'\s)ip denotes the solution of the problem 

Dtv{t) = Lris)v{t) + f{s), t>0, v{Q) = ip, (5.3) 

with 

Lr{s) := a'J{s, x)Dij + a* (s, x) A + ^^(s, x), 

for r = 1,2, Notice that the coefficients of the operator Lr{s) 

and f{s) are "frozen" at time s, thus (|5.3|) is a Cauchy problem with 
time independent data. 

We extend Theorem 12.11 as follows. 

Theorem 5.1. Let m > and k > be any integers. Let Assumptions 
\2.1\ and \2.^ hold with I > 4 + m + 4k. Let the splitting-up approximation 
Un he defined by 115.1]) or by 115. Then there exist functions u^^'^ G 
C^([0, T], W7), J = 1, 2, k, I^^ G C^([0, T], W^), such that 

Un{t, x) = u{t, x) + 6u^^^ {t, x) 

+5\'-^\t,x) + ... + S''u^''\t,x) + R^^\t,x) (5.4) 

for all t G Tn, x G M'^, anc? n > 1. T/ie functions u(^\ J = 1,2,.. .,k, 
are independent of n, and 

sup ||i?i')(t)|U,p<iV<5'+^ 

iG[0,T] 

/or a// n, where N depends only on k, d, di, K, m,p, T. 

Clearly Theorem 15.11 implies that we can again accelerate the con- 
vergence of the sphtting-up approximations by considering 

k 

Vn{t, ■) := ^ bjU2Jnit, ■), tE T„, 
3=0 

where U23n for all j = 0, 1, .., are defined by either ()5.H) or by ()5.2p . 

Theorem 5.2. Let m > and k > be any integers. Let Assumptions 
lO and\KE hold with I > 4 + m + Ak . Then for all n > 1 

max \\vnit) - u{t)\\m,p < iV5^'+\ 
where N is a constant, depending only on k,d,di, K,m,p,T . 

Hence by Sobolev's Theorem on embedding into we imme- 
diately get the following result. 
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Theorem 5.3. Let m > and k > be any integers. Let Assumptions 
\2.1\ and \2.^A hold with I > 4 + m + 4A; . Let s > be an integer such 
that m > s + d/p. Then for all n > 1 

max sup 

where N is a constant, depending only on k,d,di, K,m, s,p,T . 

We prove Theorem 15.11 bv adapting the proof of Theorem 12. II to the 
time dependent case. If u„ is defined by 1)5.11) . then we consider the 
problems 

dvit) = {L{Aoit))vit) + /(Ao(t))) dAoit), v{0) = Mo, (5.5) 

di 

dw{t) = {Lr{AQ{t))w{t) + fr{Ao{t))) dAr{t), w{0) = Uq, (5.6) 

r=l 

where Ao{t), Ai{t), A2{t), Ad^{t) are defined by (Q and (jOD, and 
L{Ao(t)), Lr{Ao{t)) mean that we substitute Ao(t) in place of the time 
variable t of the coefficients of L, L^. 

If Un is defined by ()5.2|) then we consider problems ()5.5|) and ()5.6|) 
with absolutely continuous functions Aq, Ai,...,^^^, defined by the fol- 
lowing requirements: 

Af,(0) = 0, At is periodic with period {di + 1)5, 

Ar{t) = lir,r+i]{t/5), t G [0, (rfi+l)(5] for r = 0, 1, . . . , di . (5.7) 

By virtue of Theorem 13.11 equations (|5.5|) and ()5.6p admit unique 
solutions V and w, respectively. Clearly v,w E C^([0, d T], Wp), and 

v{d't) = u(t), w{d't) = Unit) for all t E Tn, 

where d' = di if Aq, Ai, A^^ are defined by ()4.2p and ()4.3p . and 
d' = di + 1 if Aq, Ai, ...,Aii^ are defined by ()5.7j) . Therefore, our aim 
is to get an equality like ()5.4j) with v and w in place of u and Un, 
respectively. 

We treat the cases of two approximations simultaneously and warn 
the reader that, in order not to repeat the same arguments twice, we 
are going to use the same notation for some objects that have different 
meaning in each case. From now on d' denotes di if we consider the 
splitting-up approximations m„ defined by (jSHj), and it denotes di + 1 
in the case of Un defined by ()5.2|) . We keep the notation Af for the set 
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of all multi-numbers a = a\a2---oij for G {0, 1, 2, d\\ and integers 
j > 1. We use also the numbers Cq, and the functions i?Q,, defined by 
(031), and (gSl), with rf' in place of dx in Observe that 

as is easy to see Lemma 14.11 still holds with d! in place of d\ in its 
formulation. 

Let 7?./ and "Kf denote the solutions of the problems 

duif) = {Lu{t) + f{t)) dt, m(0) = 0, 

and 

dv{t) = {Lv{t) + f{t)) dAo{t), v{0) = 0, 

respectively, where L := L{Ao{t)), the operator obtained from L by 
the substitution of Ao{t) in place of t in the coefficients of L. Notice 
that TZ depends on n when Aq is defined by ()5.7|) . Notice also that 

nf{t, ■) := (7^/)(ylo(^), ■) = nf{t, ■), (5.8) 

where f{t, ■) = f{Ao{t), ■). Let Qa/ denote the solution of the problem 

dv{t) = Lv{t) dAo{t) + f{t) dB^{t), v{0) = 0. 

We modify the definition of L^, fa, used for time independent oper- 
ators and free term, as follows: Lq := 0, /o := 0, 

for 7 = r G {1, 2, rfi}, and 



fyO '■— Lfy /-),, f^r '■— Lyfr (5-9) 

for r = l,2,...,di, where := Dff^, and denotes the differential 
operator which we obtain from by taking the derivative in t of its 
coefficients. As is easy to see, and are well defined if2(|7|— 1) < /. 
We use the notation and f-y for the operator which we obtain from 
L-y by substituting Ao(t) in place of t in its coefficients, and for the 
function obtained from by the same substitution, respectively. Then 
we have the following counterpart of Lemma 14.21 

Lemma 5.4. Take some functions 

h G Lp{[0,d'T],Wl), K e Lp{[0,d'T],Wi), r = 0, 1, ...,d,. 
Let u be a solution of the "equation" 

di 

du = hr dAr = hr dAr 

r=0 
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With n(0) e Wp. Assume that Lu E Lp{[0, d'T],Wp) . Then for any 
aeAf 

TZ{Bah) = CaoTZh + 6Qaoh, 

QaU = TZ^CaoLu — CarK) + SQaoLu — SQarK + BaU. 

The proof of this lemma is an obvious modification of that of Lemma l4.2l 
Next, let us use the notation 

wp = Lpw + fp. 

bmce w e C^([0, d'T], ly^), the functions Wj3 are well defined for 2\f3\ < 
/. Under the same condition the coefficients of Lp and fp have the first 
derivative in time and these derivatives are under control. Furthermore, 
dwi3 = hr dAj., where, as long as 2|/3| + 3 < /, the functions 

ho = {LLp - Lpq)w + Lfp - /^o, 

hr = L{LrW + fr), r = l,...,di, 

are bounded -valued functions on [0,d'T]. 

Then in the same way as Lemma 14.31 Proposition 14.41 and Lemma 
14.51 are obtained by the aid of Lemma 14.21 using Lemma 15.41 we get 
their counterparts, formulated as follows. 

Lemma 5.5. Let a,l3eN. If2\(3\ + 3 < / then 

Proposition 5.6. Let k, > be an integer and / > 2k + 3. Then 

i=l \a\=i i=l |a|=i+l |a|=K+l 

(5.10) 

Lemma 5.7. If k > is an integer and a, P E Af and 2{\P\ + k) + 1 < I, 
then 

K 

n{B^Wp) = 2Z Ca07^^^/37 

i=0 |7|=i 

K 

i=l |7|=i— 1 |7|='* 

where for any multi-numbers fi, v 

c^-yTZw^^ := CyTZw^, ^ B^^w^^ := B^w^, 

|7|=0 |7|=0 
|7|=0 
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In order to iterate equation we introduce the following class 

of indices. We say that 

P = Y ■■=ll'Y2'-l7 (5.11) 

is a graded multi-number of length \j3\ := j + Ui + 1/2 + ■■■ + uj, ii 
■ji G {1, 2, di}, z/j > is any integer for i = 1,2, where j > 1 is 
any integer. If i/j = for some i, then we also write 7^ in place of 7^^^ 
in 1)5.111) . Let /C denote the set of all graded mult i- numbers. For each 
(3 = Y = 7r72'--7? e /C of length \(3\ < I + 1/2 we introduce the 
following operators and functions: 

L, = W :=(-l)'"'"'4?^^^^---4?^ 



fp = := (-1)1^1-147) ■ ... ■ L^:f>f^^\ (5.12) 

(s) (s) 

where fr '■= D'^fr, and Lr denotes the operator which we obtain 
from Lr by applying the derivation Df to each of its coefficients. By 

definition fi^^ = fr and = L,.. It is easy to see that for /3 G A/", 
when /? = G /C, definitions ()5.12|) are consistent with ()5.9|) . 

Lemma 5.8. The following statements hold. 

(i) Let P,-f eK. be such that \p\ + I7I < 1 + 1/2. Then 

LpL^ = —Lp-y, Lpf^ = —fis-y- 

(a) Let a E M be such that p := |a| < 1 + //2. Then there exist 
constants 0(7) = c(a,7) G {0,±1,±2,..} defined for a// 7 G /C with 
I7I = p, such that 

La= c(7)L^, /„ = c(7)/^. (5.13) 

7e/C,|7|=p 7e/C,|7|=p 

Proof. Part (i) follows immediately from the definition ()5.12|1 of Lp, ff^. 
Part (ii) clearly holds for a = and a = r G {l,...,di}. Assume that 
equations ()5.13p hold for some a G A/", |a| < 1 + //2. Then 

Lar = -LaLr = - ^ c{(3)LpLr = ^ c{(3)Lpr, 

/3e/C,|/3| = |a| ()&K.,\P\ = \a\ 

far = -Lafr = - ^ c{[3)Lpfr = ^ c{P)f(3r 

/3G/C,|/3| = |q| l3eK.,\l3\ = \a\ 

for r G {1,2..., di}, and 

di 

LaO = LLa - La = ^ ^ c{P)LrLf3 - ^ c{Y)L^-^, 
r=l l3£M,\l3\=p Y&M,\Y\=P 
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di 

r=l l3&M,\P\=p Y&M,\Y\=P 

Hence by using assertion (i) and noticing that 

1/^1=1 ImI=i 

we get equations (j5.13|) for ar with r = 0, 1, (ii. Thus the induction 
on the length of a completes the proof. □ 

For (3 & IC we write and //j, when the time change Ao{t) is 
done in the coefficients of Lp and in fp. We set := L-yW + f-y 
for 7 G /C, I7I < 1 + //2. Notice that Lemma 15.81 has an obvious 
translation in terms of these functions. Namely, by Lemma 15.81 (ii) 
for every a G A/" such that p := \a\ < 1 + I /2 there exist constants 
0(7) = c(a, 7) G {0, ±1, ±2, ..} defined for all 7 G /C with I7I = p, such 
that 

7G/C,|7|=p 

For every integer z > 1 we introduce finite sequences a := /?2, A) 
of graded multi-numbers /5j G /C, and we set |cr| := + |/32| + ... + 
We also introduce the empty sequence e of length |e| := 0. The set 
of all these sequences is denoted by X. For a = {Pi, P2, Pi) with 
|cr| < 1 + //2 we define 

So- '■= T^Lp^TZLp^ ■ ... ■ TZLp., '■= TZLpJZL^^ ■ ... ■ TZL^-, 

and for \a\ = we set 

■ —— "l^f . —— 7^. 

Notice that for any g G Lp([0,T], W?'") 

S„g{t,-) = {S^g){Ao{t),-), (5.14) 

where g{t, ■) := g{AQ{t), ■). This follows from (j5.8p by induction on \a\. 
In order to formulate the counterparts of Lemma (4.71 and Proposition 
14.81 we use the following sets 

A{i) = {(or,/?) : a G J,/3 G /C, |a| + \p\ < i}, 
B{z,j) = {{a, /5) : a G AT, /3 G /C, |a| = z, < j}, 

i 

B*{z,j) = [j B{zi,j). 
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Remember that if is a collection of functions indexed by a parameter 

* 

u taking values in a set A, then ^ g,^ means any linear combination 
of g^, with coefficient independent of the argument of g^ and of S. 

Lemma 5.9. Let a E 2 , ^ > be integers, and a G A/",/? G fC. 
Assume that 

2{\a\ + + + /i + 2 < /. 

Then 

i=0 A(\a\ + \l3\+i) 

K 

i=l B{\a\+i,\a\ + \f3\+i~l) 

Proof. We can derive this lemma from Lemma 15.71 in the same way as 
Lemma 14.71 is proved. We need only use the sets K. and X in place of 
A4 and J', and the operators TZ, L^, S„ for i/ G /C, cr G X, in place of 
71, Ly So-, ioT u E J^, (J E J , respectively. □ 

Proposition 5.10. Let k,m>0 be integers, and 

A + m + Ak<l. 

Then for any j = 0, 1, k we have 

j k 

i=l A{2i) i=j+l A{i+j+l) 

k 

i=l B*{i,i+j) 

where Vp := Lpv + fp. 

Proof. The proof of this proposition is a straightforward translation of 
the proof of the corresponding proposition, Proposition 14.81 in the time 
independent case. To make this translation we use the sets K, and X in 
place of and J , and the operators 1Z, L^, for G /C, cr G X, in 
place of TZ, Ly S^, ioi v E M., a E J , respectively. □ 

Now we can finish the proof of Theorem 15. II as follows. Taking j = k 
in Proposition 15. lUl we get 

j 

w = V + ^5' ^ c{a,P)SoVf3 + ^ c{a,P,5)Ba^Wf3^+rs, (5.15) 

i=l A{2i) B*{k,2k) 
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where c(cr, c(a,/5,5) are certain constants, and is a function in 
C^([0,T], ly™) for each 6, and 

sup sup 6~^'+'^\\rs{t,-)\\m,p<N. 

telO,d'T] n,S=T/n 

Observe that in contrast with ()4.21|) the functions v and ScrVf3 in ()5.15|) 
may depend on 6. To proceed further, define Rn\t,x) := rs{d't,x), 
and 

:= ^ c(cr,/3)S'^M/3, i = l,2,..,k, 

A(2i) 

where := Ljsu + fp. Then by virtue of equahty (|5.14|) and the fact 
that v{t) = n(y4o(t)) from equation (jS.lSj) we get 

wit,-) = uiAoit),-) + Y,S'^^'\Mt),-) 

4 = 1 

+ c{a,P,6)B^,{t)wp,{t,-) + R^^\t/d',-). 

B*(k,2k) 

Substituting here d't in place of t we get the required representation 
()5.4j) by taking into account that 

w{d't) = Unit), Ao{d't) = t, Ba,{d't) = ^teTn. 

Remark 5.1. Let 1 < J < c?i. Then Theorems 15.11 and 15.31 hold 
also when the operator S^''^(t"_|_;^) is replaced with 'B^p(t'^) for every r = 
l,2,...,j in the definition ()5.2j) of the splitting-up approximation 
To see this we need only repeat the proof of the previous theorem with 
Aj in place of in equation ()5.6|) and with Aq and Aj interchanged 
in (|K3|) . 
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